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Abstract 



We use the Liouville-von Neumann (LvN) approach to study the dynamics and 
the adiabaticity of a time-dependent driven anharmonic oscillator as an example of 
non-equilibrium quantum dynamics. We show that the adiabaticity is minimally 
broken in the sense that a gaussian wave packet at the past infinity evolves to 
coherent states, however slowly the potential changes, its coherence factor is order 
' of the coupling. We also show that the dynamics are governed by an equation 

of motion similar to the Kepler motion which satisfies the angular momentum 
Qh! conservation. 
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1 Introduction 



The quantization of a time-dependent harmonic oscillator has been important in many 
branches of physics such as the coherent states and the squeezed states the quantum 
statistical properties of radiation |2| , the information theory || , and the quantum theory 
out of thermal equilibrium |3J . As a more realistic system anharmonic oscillator has been 
widely studied for various purposes. In this direction, considerable efforts have been 
devoted to constructing systematic procedures to improve the variational approximation 
method of quantum mechanics and quantum field theory |, 0]. These methods 
give appropriate effective potentials which are adequate for describing time-independent 
systems. However, if one is interested in the explicit time evolution of time-dependent 
systems, above methods are limited in their uses. 

The path integral approach || was developed to study several kinds of time-dependent 
oscillators with perturbative external forces. The LvN approach to the generalized 
driven oscillator |§ was developed and its geometric phases and wave functions were 
discussed JK| . The LvN approach to the time- dependent harmonic [|ll[] and anharmonic 
oscillators [|l^] have also been developed. A merit of LvN approach is that it can be di- 
rectly applied to both time-independent and time- dependent quantum systems without 
any modification. 

In the case of a driven harmonic oscillator, it was shown 1THI that, the transition 



amplitudes between different states are zero in the adiabatic limit of a driving force of 
the form Fe~ l l T . The adiabatic assumption, that this phenomena holds for a general 
type of quasi-static potentials is usually accepted but no concrete proof has been given 
in general situation. In this paper we test this assumption by checking the evolution of a 
gaussian ground state of the time-dependent driven anharmonic oscillator. We formulate 
the quantization of a time-dependent driven anharmonic oscillator with general time- 
dependent mass, frequency, and driving force. The present formulation can also be 
applicable to the study of phase transition using the inverted mass square [[[^ . 

The organization of this paper is as follows. In Sec. II, the LvN approach is elabo- 
rated to be applicable to the driven anharmonic oscillator and we find the Fock space 
and expectation values of some operators. Sec. Ill is devoted to study driven anhar- 
monic oscillator system which is appropriate to analyze the adiabatic assumption. We 
show that the initial gaussian states evolve to final coherent states with its squeezing 
varies in time by 0(A). In Sec. IV, we consider the weak coupling limit of the oscillator 
in Sec. Ill and we show that the squeezing becomes time-independent in the adiabatic 
limit. Therefore, the adiabatic assumption is broken minimally. In Sec. V, we summa- 
rize our results and give some discussions. We have two appendices which present some 
calculational details. 
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2 Time-dependent driven anharmonic oscillator 



The quantization of a time-dependent anharmonic oscillator using the LvN approach 
was tried in Ref ||12||. It was shown that their approach is equivalent to the variational 



perturbation method |T4| and the mean-field, Hartree-Fock method. In this section, we 
extend the LvN method to be applicable to the driven anharmonic oscillator, whose 
hamiltonian is given by 



H(t) 



m 2 {t)X{t) 



2 + -m(t)u 2 (t)q 2 (t) - ^m(t)F(t)q{t) + - q\t), 



(1) 



where we have introduced the mass m(t) in front of the coupling and the external force 
for simplicity of the later equations. The mass dimensions of / and A are 3/2 and 3, 
respectively. 

The LvN approach resorts to find a pair of invariant operators 



A = i[v*(t)p(t)-m(t)v*(t)q(t)+x*(t)}, 
A ] = -i [v(t)p(t) - m(t)v(t)q(t) + x(t)) , 

which satisfy the LvN equation, 



(2) 



(3) 



ait 



dt dt 
and the standard commutation relation 



-t[A,H} = 0, 



[i,it] = l. 



(4) 



The commutation relation identifies A and A> as annihilation and creation operators 
respectively and it gives a Wronskian condition for v(t): 



h[v(t)v(t)-v(t)v*(t)]=i, 



(5) 



where v(t) = v(t)/ \Jm(t). Later in this paper, we use the unit which makes h = 1. 

The ground state is the one annihilated by A and the n th number state is obtained by 
applying A' n times to the vacuum state: 



A|0> 
\n) 



0. 



[it]r 



|0). 



(6) 
(7) 



The wave function of the ground state to this order of approximation is given by 

1/4 



i 



2nv*(t)v(t) 



cxp 



m(t)v*(t) , x*(t) 

i : , y q 2 - 1— y-q 



2v*(t) 



v*(t) 



2 



The LvN equation reduces to the following operator equations: 
x{t) = -F{t)v{t), 



with the new frequency u> 2 (t) = uJ 2 (t) 



3! 

d 2 ^m(t) 



(9) 
(10) 



Differentiating Eq. ( |TUD with 



y/m{t) dt 

respect to q(t) and taking the vacuum expectation value, which is equivalent to the 
variational approximation [12|, one gets 

m(t)X(t) 



0(f) 



<0|g 2 (f)|0> 



v{t) = 0. 



From the first basic comparison theorem of ordinary differential equation, one can natu- 
rally conclude that the solution for v(t) should be oscillatory if Co 2 and m W A W (0|g 2 (t)|0) 
is positive definite. The quantum operators q(t) and p(t) can be written in terms of the 
creation and the annihilation operators as 



q(t) = v(t)A + v*(t)A* + 2%[v(t)x*(t)}, 

p{t) = m(t){v(t)A + v*(t)A ] + 2Q[v{t)x*{t)]} 

where ^ means the imaginary part and 



(12) 



x(t) 



F(t)v(t)dt. 



The vacuum expectation values of operators are then given by 

(0|g(t)|0> = 2Z[v(t)x*(t)\, 
(0\p(t)\0) = 2m(t)Z[v(t)x*(t)}, 
(0\q 2 (t)\0) = \v(t)\ 2 + 4\%[v(t)x*( 



(13) 



(14) 



(0\p 2 {t)\0) = m 2 (t)\v(t)\ 2 + Am 2 (t) \%[v{t)x*{t)]\ 2 . 

The quantum mechanical dynamics of the forced anharmonic oscillator system are then 
contained in one classical equation of motion ([TTJ) , which becomes 

i(t)+u 2 eff (t)v(t) 

where 

A(t) 







(15) 



\v(t)\ 2 + \v(t)x*(t) -v*(t)x(t)\' 



(16) 



In summary, we have reduced the problem of quantization of driven anharmonic oscillator 
into the problem of solving one classical differential equation flT5"D. All information of the 
dynamics are encoded in a single function v(t). Once the explicit form of v(t) is known, 
we have time evolution of the quantum operators pit) and q(t), and time evolution of 
any other physical quantities are derived from these operators and the states (^). 
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3 Anharmonic oscillator with a linearly growing ex- 
ternal force 



In this section, we present the evolution of the driven anharmonic oscillators with con- 
stant and linearly increasing forces. Let us first analyze the case of time-independent 
driven oscillator. We assume ui(t) and F(t) to be time- independent. The mass m(t) 
and the frequency u(t) may be dependent on time, whose time dependency leads to no 
trouble as long as they are determined to give time-independent Co. 
An exact solution of Eq. (pl|) can be found by assuming 



x(t) = —Frexp(-iQt), Vo(t) = — iflr exp(— iQt), 



;i7) 



where F and r are constants. The value of r should be given from the Wronskian 
formula (H), and the resulting solution is 



v (t) 



1 



-iQt 



2Q 



Substituting flT|) to flU) one gets the gap equation 



n 2 



~ 2 A AF 2 
U + 4fi + 2^ 



(19) 



m{t)Q 2 
m(t)F 



which has only one positive root for Q. Note that this gap equation exactly reproduces 
the result of Ref. [12] if one turns off the external force (F = 0). 

The expectation values of operators with respect to the ground states defined by (Bp are 



(0|g(t)|0> 

(0|p(t)|0> 

(0|g 2 (t)|0) 
(0\p 2 (t)\0) 
(0|(Ag) 2 (Ap) 2 |0) 



20? 
1 



2m(t)n 
m(t)tt 



m{t) 



+ 



1 + 



m(t)Q 4 ' 
m 2 (t) 



(20) 



1 + 



4H 2 m 2 (t) 
m 2 (t) 



+ 



m 2 (t)F 2 
4Q 4 m(t) ' 



4n 2 m 2 (t) 



Note that the state |0) satisfies the minimal uncertainty relation. In fact, |0) is the 
minimum energy state of the Fock space. The force dependence on the uncertainty 
enter only through the frequency indirectly. 

Now let us go into the main problem in this section: a time-dependent anharmonic 
oscillator with the driving force linearly increasing to a maximum value /. For simplicity 
of the problem, we remove the time-dependence of the parameters except for the force. 
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As a typical type of linear driving force, we consider the force of the form, 



F(t) 



For notational simplicity, we use t. 
in Eq. ([13]) then becomes 



0, t < , 
ft 

V, < t < A , 
A 



(21) 



t > A. 



0) and t + (= A) when it is convenient. The x(t) 



x(t) 



F(t')v(t')dt' 



0. 



t < , 

±(w*e+i£rt£dem, o<t<A, (22) 

-/ f v(t')dt' + [ dt' t dt"v(t"), t > A. 
Jo A Jo J o 



Now let us solve Eq. ([15]) at each regions of time. For t < 0, the external force 



vanishes, and we find an exact solution of Eq. (|T5|) which describes the evolution of the 
ground state of the oscillator: 



v-(t) 



1 



-iQ-t 



(23) 



where f2_ denotes the value of Q(t) in Eq. ( P^|) with no external force, F — 0. 

For t > 0, we cannot find exact solution of the full equation Eq. (0). Instead, we try 
to find approximate solution based on the WKB-type solution. x(t) is proportional to / 



and the 



vx 



V x\ 



term in Eq. ([16]) would be 0(\f 2 / 



Therefore, we introduce 



the approximation for f 2 \ Q 6 and find the solution for v(t) up to linear order in 
Xf 2 /Q G . Any types of forces can be treated approximately if it satisfies Xf 2 8u 6 . Note 
that since it restricts only the product of the coupling and the driving force squared, 
the driving force need not be small. The present method used in this paper can cover 
any time-dependent external forces which approaches to a finite value at t — > oo if 
the coupling is small enough, which is usually required even in the absence of driving 
force. One may try to find an approximate solution of Eq. fll5|) by considering linearized 
equation for v(t) — V- (t) and by using (t) as a zeroth order solution. But this approach 
based on v_ (t) fails to give appropriate description of the system because the size of the 
perturbation linearly increases as QA at t = A, which gives divergent result in the 
adiabatic limit (A — > oo). In fact, more appropriate zeroth order solution is the WKB 
type: 



v (t) 



-ii>(t) 



2Q(t) 



(24) 



where ijj(t) = J tt(t')dt', and the frequency Q(t) is determined by the gap equation (jT 
in which the driving force F(t) is now time-dependent. Let us set the two values of 
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frequencies, Q(0) and Q(A) as Q T , which are determined by 

Q 2 _ — u) 2 + — ^— , ^+=^ + ^- + ^. (25) 



For later uses let us define two dimensionless constants: 

A \f 2 



(26) 



16fii' 8fi 6 (l + 2z/)' 
The phase in the intermediate time between < t < A can be calculated up to 0(e) as: 

ip(t) = / n(t')dt' = Q^t + 4r> < t < A. (27) 

Jo 3 A J 

Now we replace v(t) in \v(t)x*(t) — v*(t)x(t)\ 2 term of Eq. C[X6f) by its zeroth order 
solution v (t). Note that we do not replace |"5(£)| 2 term by its zeroth order solution, 
since we now expand the solution in e = Xf 2 /(8Q S ). Therefore |t>(t)| 2 term is zeroth 
order in this parameter. We calculate all physical quantities up to this order. The 
equation of motion fll5|) then becomes 

fi(t) + u> 2 (t)v(t) + ^\v(t)\ 2 + j^P®}*® = °> ( 28 ) 

where f(t) is a dimensionless function defined by 

fit) = -i^±[v(t)x*(t)-v*(t)x(t)\ 
Q 2 (t) 

- -i^[m^(t)-v*{t)x(t)] Ht)=Vo(t) (29) 

r o, t < o , 

' [ifj(t) -sin^(t)] , 0<t<A, 



Afi(t) 



sin(fi+A/2) „ / A\ 

y cos 0+ ( t - — , t > A. 



fi+A/2 + V 2 

During the calculation, one may find the following term 
sin [(V»(A) - Of A)/2] 



Afi/2 



cos[n+t + n+A/2-V(A)/2], (30) 

added to the last line of Eq. fl2"9"|) . This term can be removed by the following argu- 
ments: First, consider the case Af2_ C 1/e. Then, the argument of the sine function, 

— — — = — efi_A in Eq. fl5C|) is 0(e), which makes above terms to be 0(e). On 

2 3 

the other hand, if A 3> 1, the denominator in Eq. ( p0[ ) diverges. Therefore, this term 
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can be set to zero safely to 0(e°). There appears many similar terms in the calculations, 
which we simply omit without explicitly referring to it. 

Because the driving term in Eq. ( p^) is linear in v(t), the driving term in the potential 
which gives Eq. (p§D is quadratic in v(t). If there is no higher power of v(t) in the 
potential one cannot prevent the solution resonating to the driving source and becoming 
divergent, which is unwanted result. Therefore, the higher power term in the potential is 
necessarily needed, which is the A|5(t)| 2 term in Eq. (p8|) . Even though we cannot solve 
this non-linear time-dependent differential equation exactly, it is possible to obtain the 
general properties of the solution by separating v(t) into its size r(t) and phase 4>(t): 



v(t) = r{t)e 



i<t>{t) _ 

Then the equation of motion ( p8| ) is reduced to the following two equations: 

L 2 




m ~ r»(t) 
r 2 (t)0(t) = L 2 , 

where L is a constant of motion determined by 
L 



r 2 (t)0(t) = r 2 (A)0(A) = r 2 (O)0(O) = -. 



(31) 

(32) 
(33) 

(34) 



Note that Eq. ([34]) is very similar to the conservation law of the angular momentum. 
Moreover, Eq. © is very similar to the radial equation of 2-dimensional rotating par- 
ticle, except that the r dependence on Eq. (|32|) is 1/r 3 rather than 1/r 2 . The effective 
potential, which gives Eq. (|32]) , can be written as 



(35) 



V±(r) + 



2fi 4 



P(t)-f z 



Li 
2r 2 



+ |r 4 + 



± 



Co 2 + 



Mlf, 



are the adiabatic 



where /+ = 1, /_ = 0, and V T (r) 
(A — > oo) potential at t < 0, and t > A respectively. If the system is quasi-static 
(QA ^> 1), the function f(t) in the potential (|35| ) monotonically increases from (t = 0) 
to 1 (t — A) in time, which define V±(r). The two stable circular orbit r = r_ (at t < 0) 
and r + (at t > A) are determined by the formula V±(r) = 0, where prime denotes the 
derivative with respect to r, and are given by 



(36) 



The evolution of (r(t),0(t)), therefore, can be seen as an orbit change of a circulating 
particle by a varying frequency with somewhat odd angular momentum dependence. 
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Now let us approximate the potential up to quadratic order in e near r = r±: 

V(r,t)~V eff (r,t) = V ± (r ± ) + lv^r ± )(r-r ± r-^^J(t)(r-r ± ) : 

where we have written the time- dependent term up to the first order in r — r± since it 
is of 0(e) already and the source 

j(t) = e(t)[e(A-t)-p(t)}, (37) 

is dimensionless. 

Let the orbit difference be 

Pi-{t) = — [r(t)-r_], 0<t<A, (38) 
e _ r _ 

Pi+(t) = — [r(*)-r+], t>A, 

A p 

where e± = —r-z is a dimensionless expansion parameter, and u)i = VI' (r±). The 
2S2±u;± 

equation of motion for p(t) is given by 

p/±(t)+a>|p /± (t) = u>i J(t), (39) 

where the subscripts ± are used for t > A and < t < A, respectively. This is just the 
equation of motion of a harmonic oscillator with a driving force, and its inhomogeneous 
solution can be written explicitly using the Green's function method: 

Pi±(t) = f ou±dt'J(t')sinu±(t-t'). (40) 
Jt± 

Note that this inhomogeneous solution satisfies pi±(t±) = = pi±(t±). We present the 
explicit form of pi±(t) in Appendix A. 

With this pi±(t) in hand, it is easy to write the solution for r(t) which satisfies the 
appropriate boundary condition at t — t±: 

r(t+) 

r(t) = r± + er± Pl ±{t) + [r{t±)-r±] cos u±{t-t±) + ^^- sin a>±(t-t±) (41) 

u± 

= r± + s±(t), 

where s±(t) denotes the time-dependent parts of r(t) in the region t > A and < t < A, 
respectively. We have omitted the sign ± in e here and later since e± are the same to 
this order. Usually, the homogeneous terms are first order in the expansion coefficient 
(e here). 

The boundary condition at t — is 

r(0) = -^=, f(0)=0, 0(0) = 0, 0(O)=fi_, (42) 



where r(0) was determined by Eq. ([|). After matching these boundary conditions we 

get 



1 3A 
L = -, r_=r(0), c^=4cu 2 + — = 4^(i + 2 z,), 



and 



r(t) 



yfWT 



[l + epi-(01- 



(43) 



(44) 



The solution for r(t) in the region t > A, which matches the boundary condition at 
t = A, is given by 



r(t) = r+ + s+(t) 



(45) 



r+ + er+p /+ (t) + -== [pi-(A) + 1] cos^+(£ - A) 

e P/-( A )„:„.- 



sinc<j + (t — A), 



where u\ = V"{r + ) = 40^.(1 + 2z/). 

In the region < t < A, the phase function can be written (See Appendix B for 
details.) as 

4>(t) = if>{t) + eS<f>.{t). (46) 

where ip(t) is given by (^) and 50_ (t) is a finite oscillating function defined in Eq. (|84j ). 
With these results we obtain -y(t) in the region, < t < A, 



i?(t) = v _(t) + ev_(t)w-(t), 



where 



«-(*) 



2f2(t) 



and u>(t) is a function of 0(1). Then x(t) in Eq. (|i~3|) is given by 

X(t): 



t 2 \ fvJt) , . 



(47) 
(48) 

(49) 



where Xi(t) is defined in Eq. (|88D in Appendix B. Now we compute the expectation value 
of operators: 

(0|g(*)|0) = ^ {[^ - sin^(t)] {1 + eto[w(t)]} - e [1 - cos^(t)] <*[«;(*)] - eO^)} , 



(0|p(t)|0> 



m 



An 2 



[1 - cos^(t)] 



, 2et 2 „ / , , .w(t)' 
1 + — + u>(t) + «- 



A 2 



(50) 



+e[fit-sinV(t)]3f \w(t) +i 



In the region t > A we get quite good approximation of the phase up to 0(e): 

t - £1 



4>{t) 



2f 2 



+ a50+(t), 



where 5(j){t) is an oscillating function of 0(1), and r = r 
Eq. (|9l|) . Then v(t) and x(t) can be written to this order, 

v(t) ~ [1 + ew+(t)] u+(t) = [1 + ew + (t)] fe- i(t -* l)/(2f2) 

sin[A/(4f 2 



i _i_ sin 2 QA/2 
1 + 6 (HA/2) 2 



z(t) 



2tf 2 fv + (t) 



-e 4r 2 — 



+ ez{t) 



(51) 
is given in 

(52) 
(53) 



A/(4f 2 ) 

where ti, w + (t), y, and are given in Appendix B. Using these we calculate the 
expectation values of operators p(t) and q(t) at t > A which shows the effects of the 
driving force. The expectation value of q(t) and p(t) in this order are given by 

(0|$(t)|0) (54) 



= 4/f 4 | 
(0|p(t)|0> 



sinA/(4f 2 ) t-A/2 

■ cos - e 



A/(4f 2 ) 



2f 



e a? 2 



y* + e i i^y + -z*(t)-z(t) 



m 



2/f 2 



1 



+e3 (w+(t) + 2if 2 w + {t) 
+ 2e/f 2 



3?(w+(t) + 2if 2 w+(t)) 
sin A/(4f 2 ^ 



sinA/(4f 2 ) . t-A/2 



; t— A _ „■ t — A 

l/e "J? 2 " —ye 



A/(4f 2 ) 

The fluctuations of g(t) and are given by 
(0|[Ag(t)] 2 |0) - ^ 
(0\[Ap(t)} 2 \0) 



cos 



A/(4f 2 ) 
t- A/2 X 
2f 2 



sin 



2f 2 



(55) 



^)| 2 = -[l + 2e^(t)] 



(56) 



™ l*(*)r=4^[l-2e«K*)l 



where is a purely oscillating function of time of order 0(1) defined in Appendix 
B (PP|). Therefore, the uncertainty to order e becomes 

1 



((Ap) 2 (Aq) 2 ) 



4 



+ 0(€ t 



(57) 



one can naturally identify that l/(2r 2 ) plays the 



If one compares Eq. ( pq) with Eq. 
role of frequency of the oscillator. 

The main effects of the driving force to the dispersions of q(t) and p(t) come in two 
fold: First, driving force makes the dispersions fluctuating [s(t)]. Secondly, it perma- 
nently squeezes the wave-packet by distorting its frequency, [l/(2f 2 )]. Note that the 
adiabatic limit, f = r + , is the least squeezing state of the system. An interesting point 
is that there are other systems with A/(4f 2 ) = nit which reveals the same least squeez- 
ing as the adiabatic limit. The interpretation on this phenomenon is given in the last 
section. 
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4 Weak coupling limit 

In this section we consider the weak coupling limit [y = A/(16fi 3 ) <C 1) of the time- 
dependent anharmonic oscillator with driving force given by Eq. (PI|) . One should note 
that there are terms of 0(1/ v) given in pi±(t) in the Appendix A, if the coupling is very 
small. The inhomogeneous solutions are approximated in the weak coupling limit as 



Pi-{t) 



Pi+(t) 



1 



4uA 2 Ql 



[cosu;_t — cos2^>(i)] 



A 2 n 2 _ 



sin 2 (fi+A/2) . (u+ + 2tt + )t - u + A . (u+ — 2Q + )(t - A) 



2*/(n+A/2) 2 
sin 2 (fi+A/2) 
' 2{Q + A/2) 2 



sin 



sin 



(5? 



(59) 



+ g"(t), 



where g'(t) and g"(t) are functions, whose integrals over time from to t are of 0(e°). 
We ignore 0(A) terms in v(t) and discuss these terms later. To this order, v(t) is given 
by 



v(t) 



v-{t), 0<t<A, 
v+{t), t>A, 



where 



«-(*) 
v + (t) 



u_(t) |1 + i(3 sin 
zs_(A) 



UJ-t 



3 i[^(i)+0_t/2] 



1 + 



+ 



MA) 



is_(A) 



,t((3 + -A)t 



(60) 

(61) 
(62) 



+ 2^+(t)e^ +{t) sin^-(t), 
where the functions if)±(t), v + (t), and the coefficients are defined by 



n + yt-A) + n + A, ^(t) = (^±-n+)(t-A). 



(63) 



v + (t) 



/ 2 sin 2 (ft+A/2) 



2n%(n+A/2) 2 

The initial condition (|45|) at t = A is given by 



/ 2 - 1 
/? ~ AW n ~ 2f2' 



S _(A) 
s_(A) 



_/3 < sin 



u;_A 



+ ^(A) 



sm 



uj_A 



-V'(A) 



-/3f2_r_ < cos 



lD_A 



+ V(A) 



sm 



a;_A 



^(A) 



(64) 
(65) 



Note that these terms can be removed by the following arguments: First, consider the 
case fiA <C 1/V. Then, the argument of the sine function, — if) (A) = f2_A (u — yj 
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in Eq. (|64| ) is 0(e) which makes above terms to be 0(e). On the other hand, if A ^> 1, 
/3 goes to zero in this case. Therefore, the homogeneous terms in the region t > A can 
be set to zero safely to this order. There appears many similar terms in the calculations, 
which we simply omit without explicitly referring to it. 

What we want to know is the expectation values of operators at t > A. Therefore 
we should calculate x(t) in this region of time. At t — A, by integrating —F(t)v-(t) in 
Eq. (|6l|) , we have 



x(A) 

In general, for t > A, 



sin^(A)/2 
Q_A/2 ' 



,#(A)/2 



z(t) = ar(A) - / / v(t')dt' = x + (t) + x(A) - x+{A), 



where 



x + (t) 



F(t)v + (t')dt' 



(66) 



(67) 



(68) 



ifv+(t) 



sin(fi+A/2) 
Q+A/2 ' 



,ifl(t- A/2) 



- 2iasmip_(t)e i ' p+{t) } . 
At t = A the following equality holds: 

iL(A) =v-(A) =v + (A) + 0(e). 

We therefore have, 

x(t) = x+ (t) + 0(e). 
The expectation value of q(t) is then given by 

/ 



(69) 



(70) 



(0|?(*)|0> 



[1 + 2a cos ip + (t) sin vp_ (t)] 



sin(fi + A/2) 
fi+A/2 



cos[fi(t - A/2)] + 2acos^+(t) sin tp_(t) 



2af 



sin ip+(t) sin ?/>_(£) 



-^-.^-A/*)] (71) 



-2asin^+(t) sin^(t)] +0(e). 

The solution is composed of the multiplication of two parts: the rapidly oscillating part, 
ip+(t), and slowly oscillating part, — ~~ A). The Eq. (|7lD is exact answer 

to all order in a for a ^> e. The correction to Eq. (|7T| ) begins only at the 0(e). It 
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is not needed to calculate the expectation value of p(t), because it is simply given by 

<p(f)> = m«. This proves that the adiabatic theorem holds in the absence of 

0(e) correction. The 0(e) corrections are given by Eq. ([54]) with the replacement g(t) 
in Eq. (|83T) by g'(t), g(t) in Eq. fl89|) by g"(£)- The adiabatic limit of weak coupling 
expansion is written explicitly in Eq. (|75|). The dispersions are given by 

(0|[Ag(f)] 2 |0) = — ^ [l -4asin^+(t)sinV'-(t) + 4a 2 sin 2 ^!(t)l , (72) 
2m\l L J 

TYlQi r 

(0|[Ap(i)] 2 |0) = — [l + 4a sin ^+(t) sin ?/>_(£) + 4a 2 sin 2 ^! (t) 

Note that the uncertainty is 1/4 to 0(a) even though the dispersions are fluctuating in 
0(a). In the adiabatic limit, the uncertainty becomes exactly 1/4, which is the minimal 
one. 

The 0(e) calculation is not important except for the adiabatic limit and the case 
AQ + = mr. In the adiabatic case we explicitly calculate each functions again up to 0(e) 
and we get, in the region t > A, 

V(t) = _^ e -^+(*-A)-#(A) ) (73) 

v 2fi + 

x(t) = ^h+te-^H, (74) 

i Z_|_ L J 
Therefore the expectation values of operators are given by 

/ 

n 



(0\q(t)\0) = 7^{l+4ecos[fi + (t-A)]}, (75) 



<0|p(t)|0> = -^Bm[n + (t-A)]. (76) 



The dispersions are 



<0|[A«*)P|0> = ^L-, (77) 

(0|[Ap(«)] 2 |0> = (78) 

The uncertainty takes the minimal value. Therefore, in the adiabatic limit of weak 
coupling expansion, the initial ground state becomes a gaussian wave-packet, which is 
very close to the new ground state except that the center of the packet is oscillating, 
whose amplitude is of 0(e). 

5 Discussions 

We have formulated the quantization procedure of a time- dependent driven anharmonic 
oscillator using the Liouville-von Neumann approach, and then applied it to the case 
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of anharmonic oscillator with linearly varying driving force. We have found a new gap 
equation for the driven oscillator and examined whether or not the adiabatic conjecture 
hold for the anharmonic oscillator. We also have calculated the weak coupling limit of 
the system. It was shown that the adiabatic conjecture breaks down at the first order of 
the coupling and the gaussian ground state evolves into a coherent state with the 0(e) 
correction of squeezing parameter. 

The weak coupling expansion is most interesting because there appears a new expan- 
sion parameter a ^> e (|63D, which are closely related with the force gradient //A. In the 
adiabatic limit of weak coupling expansion, all dynamics related to the new parameter 
a disappears and only the very small 0(e) correction remains. The expectation values 
of operators in the adiabatic limit of the weak coupling, at t > A are given in Eq. (|75|) . 
The dispersion (|77|) implies that the wave-packet has the shape of a ground state of an 
oscillator with frequency Q + . This state of the adiabatic limit is the coherent state with 
the coherence parameter of 0(e). Explicitly, the invariant operator can be rewritten as 



A = iv*(t) lp(t) - imtl- 



m 1 



m' 



+ ^fe^ (79) 

i Z_i_ 



= e in+(t-A)+#(A) B{t) + 4e ^ c #(A) , 

where B(t) is an annihilation operator of the harmonic oscillator whose center of the 
potential is fixed at with natural frequency This observation naturally shows 

that all initial states of the anharmonic oscillators are mapped into the corresponding 
states of the new vacuum except for 0(e) displacement. 



There are infinitely many ways which gives states satisfying Eq. (|7^). If A/(4r 



nn 



in Eqs. ( |54] , |55| ) , Eq. (f75|) naturally holds and the corresponding final state is very close 
to the new ground state. These arguments hold also in the weak coupling limit of the 
system Eq. Q72|). The natural interpretation of this phenomena is the following: After 
turning on the driving force, the center of gaussian wave-packet starts to move. It 
periodically passes the bottom of the potential, and accidentally, if the driving force 
stop to increase when the center stays the bottom of the potential, the packet stop to 
move thereafter. 
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Appendix A. Source integration. 

The Green's function in Eq. ( f40D can be explicitly integrated to 



Pi-{t) 



dt' 



to. 



[ip(t') — sinip(t')] 2 smu-(t — t' 



A 2 n 2 (t) 



10) 



2A 2 Q 2 



l-4£ 2 ) (l-coscj_t) + 



sin 



16£ 2 



+ 



;i-0 2 (i+0 2 

4 



sin 



uo-t + ip(t) 



sm 



~Q-t - i/;(t) 



sm 



ip(t) sin^(t) - 2V> 2 (£) 



where £ = Q_/u_. Note that this inhomogeneous solution have a quadratically increas- 
ing term in time. There is a term of 0(1/ X) in the A — > limit because (1 — 4£ 2 ) ~ 2v 
in this limit. 

After integration of the Green's function pi+(t) is given by 



Pi+(t) 



/ dt'u+J(t') sina;(t — t') 

J A 

2sinft + A/2 



Q+A/2 

sin 2 (fi+A/2) 
2(Q + A/2) 2 



■ (a+Sl+)t-uA ■ (0-Q+)(t-A) ■ (0+Q+)(t-A) ■ 

sm v v sm 1 ±r sm v + „ sm 

Z Z _j_ z 

(g-gg+Kt-A) 



(81) 



1 - coscu(t - A) + 



sm 1 2 sm 



2X 



+ 



sin (g+gn+xt-A) gin (ip-2n + )t+iDA 



l + 2x 



where x — ^+/^+- The solution ( |8~TD vanishes if the adiabatic limit, AQ — > oo, is taken. 
It also vanishes if the time interval of the driving force satisfies A = rwr/O. Note that 
there is a term of 0(1/ X), if the coupling strength vanishes, since (1 — 2%) ~ v. 



Appendix B. Calculation of v(t) and x(t). 

The phase factor <f>(t) is given by 

, , f l dt' 

Because we usually exponentiate this phase with i, one should be cautious in approximat- 
ing this integral by separating the indefinitely growing part from the small fluctuations. 
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We calculate this integral for each regions of time separately. To have a good approxi- 
mation of the integral in ( j82|) one should separate r(t) by its frequency and fluctuating 
parts. The integration of the frequency part continuously increases in time, which gives 
the phase of v(t), and the integration of fluctuating part should converge to a small 
value of 0(e) so that its exponential can be Taylor expanded. 

In the region < t < A, following the above prescription with the explicit form of 
Pi ~{t) of (|80l) , r(i) can be written as 



r(t) 



a 2 q 2 _ 



+ eg(t) 



20(t) 



53) 



where ip(t) is given in Eq. 



ip 2 (t) 

and the integral of g(t) = pi-{t) + ^ 9 ^ 2 in Eq. fl33|) 



converges to a finite number for any value of t or A. Expanding the right hand side 
of ([53|) in e the phase integral fl32|) is approximated as 

<f>(t) ~ ip{t) - 2efi_ [ g(t')dt'. 

Jo 

With this calculation we obtain v(t) in the region < t < A: 

v(t) = u_(t) + eu_(t)i«(t), 



where 



«-(*) 
w(t) 



1 



2Q2 



A 2 fi 



+ #(i) + 2ifi_ tdt'gft') 
Jo 



Then x(t) in Eq. 

x(t) = 

where 

2Tl(t) = 

+ 



^) is given by 

fV-it) r 



AQ 2 



3e #(*) 
Q 2 



1 + zfi_t - e #(t) 



1 + e A^, 



A 



exx(t), 



2i 
3TL 



t +'2 



dt' 



irl 



dt "t_ e -m") 



v(t) 



-t / v(t')w(t')dt' + dt' dt"v(t")w(t") 



In the region t > A, r(t) can be written as: 

r(t) =f[l + eg(t)}, 



»4) 
56) 



i7) 



59) 
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where 



r(t) 

eg(t) = -y- - 1 
r 

is an oscillating function of 0(e) with vanishing time average and 

sin 2 ft + A/2~ 



(90) 



r = r + 



1 + e- 



(91) 



In the adiabatic limit, the inhomogeneous solution pi+(t) vanishes. The phase variable 
4>(t) can be exactly integrated in this case: 



0(t) = 0(A) + 



* L 



a r 2 (t') 
b 



dt' 



(92) 



sm 



u){t-tp) 



+ 



2u(rl - 6 2 ) r+ + bcos Qr ^[^Z^ 



: tan 



•2 fe2 tan M 



6 . w(t-to) , 1 
sm 1 -, 



2urf 



L0(t - t ) 



-^tan^^ 



1 + \ tan 



+ finite periodic function, 

2ri 



(93) 



where b cos Cu(t — 1 ) = [sq(A) — (r x — r )] cosu;(t — A) + S2i^l smu(t — A) and O > ti are 
constants which are needed to satisfy the boundary condition at t = A. In the case of a 
non-adiabatic change, we cannot integrate Eq. ( |52"D exactly but Taylor's expansion with 
Eq. (]89f) gives a good approximation of Eq. (p2[)up to 0(e): 



0(t) = 0(A) + 



* L 



A r 2 (t') 



(94) 



t-h Jidt'g(t" 



2 jy« 2 

where 0(A) is the phase at t = A. With these v(t) is approximated as 
u(t) = [1 + ew + (t)\ v+(t) = [1 + em+(t)] f e - i( '-' l)/(2f2) , 

where 



u+(t) 



re 



-i(t-tl)/2F 2 



w+(t) = <?(*)+ i 



(95) 

(96) 
(97) 
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Using 



x(t) = - f(t')[v + (t') + ev + (t')w + (t')]dt' (91 
Jo 

rA 



+ / dt'f(t')[v + (t') + ev + (t')w+(t') - v^t') + eu_(f>-(0], 
Jo 

one finally obtains 

x{t) = -2if 2 fv+(t) 



. sin[A/(4f 2 )] ^ t i_A . 
1 . — ^ — e if 2 — eye 2^ + ez{t) 



v2\ 



where 



sin[A/(4f 2 )] lA 

1 — — e 4r^ — 



A/(4f 



A/(4f 



D i(A-t 1 )/2f 2 -iip(A) 
(2fi+) 3 /2 f 3 



(99) 



1 sin[V;(A)/2] 



3 iV(A)/2 



^(A)/2 



,(100) 



2if 2 w+(t)z(t) 



Afi 



' u_(A)xi(A) - [ v + (t')w+(t')dt'. 
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